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surface-targeted antisecretory drugs 
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Secretory diarrheas such as cholera are a major cause of morbidity and mortality in developing countries. We previ- 
ously introduced the concept of antisecretory therapy for diarrhea using chloride channel inhibitors targeting the 
cystic fibrosis transmembrane conductance regulator channel pore on the extracellular surface of enterocytes. 
However, a concern with this strategy is that rapid fluid secretion could cause convective drug washout that would 
limit the efficacy of extracellularly targeted inhibitors. Here, we developed a convection-diffusion model of wash- 
out in an anatomically accurate three-dimensional model of human intestine comprising cylindrical ciypts and villi 
secreting fluid into a central lumen. Input parameters included initial lumen flow and inhibitor concentration, 
inhibitor dissociation constant (K^) , crypt/ villus secretion, and inhibitor diffusion. We modeled both membrane- 
impermeant and permeable inhibitors. The model predicted greatly reduced inhibitor efficacy for high crypt fluid 
secretion as occurs in cholera. We conclude that the antisecretory efficacy of an orally administered membrane- 
impermeant, surface-targeted inhibitor requires both (a) high inhibitor affinity (low nanomolar K^) to obtain 
sufficiently high luminal inhibitor concentration (>100-fold K^), and (b) sustained high luminal inhibitor concen- 
tration or slow inhibitor dissociation compared with oral administration frequency. Efficacy of a surface-targeted 
permeable inhibitor delivered from the blood requires high inhibitor permeability and blood concentration (rela- 
tive to K^). 



INTRODUCTION 

Excessive intestinal fluid secretion occurs in entero- 
toxin-mediated secretory diarrheais caused by Vibrio cholera 
and enterotoxigenic Escherichia coli (Field, 1979). The 
rate-limiting step in fluid secretion is chloride transport 
from the enterocyte cytoplasm into the intestinal lumen, 
which creates the electroosmotic force driving sodium and 
water secretion (Murek et al., 2010; Venkataisubramanian 
et al., 2010; Thiagarajah and Verkman, 2012). Cell cul- 
ture and animal models (Chao et al., 1994; Gabriel et al., 
1994; Thiagarajah et al., 2004) indicate that elevation in 
cyclic nucleotides caused by bacterial enterotoxins 
activates the CFTR, a chloride channel expressed on 
the luminal surface of enterocytes. CFTR inhibition 
is thus predicted to be of clinical benefit as antisecre- 
tory therapy in diarrheas caused by bacterial entero- 
toxins (Al-Awqati, 2002; Zhang etal., 2012). 

We identified a class of small molecules, the glycine 
hydrazides and the related malonic acid hydrazides 
(MalH) , as CFTR inhibitors that target the extracellu- 
lar-facing pore of CFTR (Muanprasat et al., 2004). An 
extracellular site-of-action was suggested by patch-clamp 
measurements showing outwardly rectifying whole-cell 
currents and rapid single-channel flicker (Muanprasat 
et al., 2004) , amd proven from CFTR inhibition by mem- 
bran e-impermeant MalH-polyethylene glycol conjugates 
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(Sonawane et al., 2006). Subsequently, multivalent 
membrane-impermeant conjugates of MalH with lec- 
tins (Sonawane et al., 2007) and polyethylene glycols 
(Sonawane et al., 2008) were synthesized with < 
100 nM for inhibition of CFTR chloride conductance. 
These membrane-impermeant, nonabsorbable con- 
jugates showed antisecretory efficacy in closed-loop 
and sucking mouse models of cholera. Membrane- 
permeable, absorbable glycine hydrazide analogues 
were also synthesized for potential therapy of polycys- 
tic kidney disease, in which cyst growth is CFTR de- 
pendent (Yang et al., 2008). 

Nonabsorbable CFTR inhibitors are potentially useful 
for antisecretory therapy because of their minimal sys- 
temic exposure. A glycine hydrazide analogue, iOWH032 
(de Hostos et al., 2011), with modest CFTR inhibition 
potency {K^ of ^8 pM) , is currendy in early-stage clinical 
trials for cholera therapy. A natural product, crofelemer, 
which is thought to act by inhibition of CFTR and cal- 
cium-activated chloride channels (Tradtrantip et al., 
2010), is also in clinical trials (Cottreau et al., 2012). Cro- 
felemer consists of large proanthocyanidin oligomers 
that are predicted to be membrane impermeable and 
thus externally acting. 
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A potentially important though largely ignored con- 
cern in the application of externally targeted inhibitors 
is washout as a result of convection caused by rapid 
fluid secretion. Convection effects cannot be addressed 
easily in animal models because of differences from 
humans in ctypt-villus geometry and fluid secretion 
rates, as well as the lack of suitable animal models of 
secretory diarrhea. As diagrammed in Fig. 1 A, the intes- 
tine contains a dense array of long, narrow cylindrical 
crypts that secrete fluid into a central lumen. Because 
of fluid convection, crypt fluid secretion is predicted to 
reduce inhibitor concentration at the enterocyte sur- 
face and hence reduce inhibitor efficacy. Further, if 
therapeutic inhibitor concentration is not sustained in 
the intestinal lumen, its antisecretory effect would di- 
minish rapidly unless inhibitor dissociation from its tar- 
get at the enterocyte surface is slow compared with the 
time between oral doses. Here, we model CFTR inhibi- 
tor washout to predict the conditions (concentration/ 
target affinity, dissociation rate) under which an extra- 
cellularly targeted antisecretory drug could be effective 
in reducing intestinal fluid secretion. 

MATERIALS AND METHODS 

Model formulation 

Fig. 1 A diagrams a cross section of a generic intestinal villus- 
crypt unit, as found in mid-jejunum, which consists of a long, 
narrow aqueous compartment bounded by a cylindrical layer of 
epithelial cells. The crypts project in an outward radial direction 



from the intestinal wall, and the villi project in an inward radial 
direction. The inhibitor target, CFTR, is heterogeneously dis- 
tributed in the lumen-facing (apical) membrane of epithelial 
cells along the crypt-villus axis, with its extracellular surface ex- 
posed to the luminal fluid. Based on published anatomical data 
in human male adult mid-jejunum (Loehry and Creamer, 1969; 
Marsh and Swift, 1969; Trbojevic-Stankovic et al., 2010), typical 
crypt length is 150 pm and inner diameter (of the aqueous lu- 
men) is 20-25 |im; typical villus length is 350 pm and inter-crypt 
spacing is ^100 pm. Volume flux across the epithelium at loca- 
tion z along the crypt-villus, Jv(z) , depends on CFTR density, the 
magnitude of the secretory stimulus, and the concentration of 
inhibitor contacting the villus-crypt wall. In cholera, typical 
fluid secretion rate in mid-jejunum is ^4.5 ml/h per cm intes- 
tine (Banwell et al., 1970; Van Loon et al., 1992; Bearcroft et al., 
1997), which, together with a crypt density of 12,500 cm"^ 
(Trbojevic-Stankovic et al., 2010) and a jejunal diameter of 4-5 
cm, corresponds to a single-crypt fluid secretion rate of ^4 x 
10"' cmVmin (~7 x 10"^ pL/cm^/s at the crypt surface). Fluid 
viscosity in the crypt lumen is ^10-fold greater than that of water, 
as measured by photobleaching and dye polarization (Naftalin 
et al., 1995) . Material properties and baseline parameters are listed 
in Table 1. For a single cylindrical crypt-villus unit, the com- 
puted parameters include inhibitor concentration Cj (r,z) and 
velocity Vj (r,z) , in which cylindrical symmetry allows specifica- 
tion of location by radial dimension r and axial position z. 

Fig. 1 B shows the generic three-dimensional geometry of the 
intestine, with crypts and villi projecting radially from the cylin- 
drical intestinal wall. The diameter of the intestinal lumen is 
orders of magnitude greater than crypt length. Fluid entering 
the intestinal lumen has a specified inhibitor concentration and 
flow velocity, which may vary in time (because of intermittent 
oral administration). Because the intestinal lumen is cylindri- 
cally symmetric and large compared with cr^ypt dimensions, dif- 
fusion-convection for the geometry shown in Fig. 1 B can 
be simplified to that of a truncated radial wedge of fluid with 




Figure 1. Model of drug convection- 
diffusion in intestinal crypt-villus units. 
(A) Schematic of epithelial cell-lined crypt- 
villus units (lengths, L^rypt and L^nus). Fluid 
secretion into the kimen produces convec- 
tive (upward) solute transport. (B) Three- 
dimensional model of radially oriented 
crypt-villus units projecting from a cylin- 
drical lumen. (C) Section near the intes- 
tinal wall showing crypt-villus units and a 
wedge of luminal aqueous-phase volume. 
Cylindrical symmetry reduces computation 
complexity. An example of pseudocolored 
inhibitor concentration profile is shown. 
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attached ciypt-villus units as diagrammed in Fig. 1 C, as described 
in Results. 



V ■ Vi = 0, 



(3) 



List of parameters 

Parameters for single-crypt/villus computations: 

D, diffusion coefficient in crypt/villus lumen (cm^/s) 

Jv(z), volume flux across the epithelium at location z along the 

crypt (pL/ cm^/ s) 

Jv°(z), volume flux across the epithelium in the absence of inhibi- 
tor (pL/cm^/s) 

Xj„ total volume flux from crypt (pL/s/ crypt) 
Ci(r,z) , inhibitor concentration at location r,z (pM) 
Vi(r,z), vectorial flow velocity at location r,z (cm/s) 
Pi(r,z), pressure at location r,z (Pa) 

AP, pressvire difference between Pi(r,z) and reference pressure in 
the intestine (Pa) 

Co, inhibitor concentration at outer boundaiy (pM) 
Jb(z), transepithelial volume flux of (permeable) inhibitor 
(mol/cmVs) 

Jb°(z), initial transepithelial volume flux of inhibitor (mol/cmVs) 

Pinh. transepithelial inhibitor permeability (cm/s) 

Cb, inhibitor concentration in blood (pM) 

K^, inhibitor-CFTR binding dissociation constant (pM) 

k.1, dissociation rate constant (niin"') 

ki, association rate constant (pM"' min"') 

p, fluid demity (kg/ m^) 

17, fluid dynamic viscosity (Pa ■ s) 

At, computation time step (s) 

L, crypt length (pm) 

d, crypt diameter (pm) 

fb(z), occupied fraction of inhibitor binding site 

Additional parameters for multi-crypt/villus computations: 

Dl, diffusion coefficient in the intestinal lumen (cm^/ s) 

Dm, diffusion coefficient in the mucosa (crypt/villus region) 

(cmVs) 

77l, fluid dynamic viscosity in the lumen (Pa ■ s) 

T)M, fluid dynamic viscosity in the mucosa (Pa • s) 

Ly, villus length (pm) 

Lc, crypt length (pm) 

dy, villus diameter (pm) 

dc, crypt diameter (pm) 

dL, intestinal lumen diameter (pm) 

U(x,y,z), inhibitor velocity in the intestinal lumen (cm/s) 

Umeani mean inhibitor velocity in the lumen (cm/ s) 

s, distance from the center of the lumen (cm) 

So, distance from the center of the lumen to the top of the villus (cm) 

Model computations 

Models of convection-diffusion in single crypt-villus units (two- 
dimensional) and the full intestine (three-dimensional) were im- 
plemented using COMSOL Multiphysics (version 3.4; COMSOL). 
The model is specified by a diffusion equation describing inhibitor 
convection-diffusion coupled with the Navier-Stokes equation de- 
scribing the fluid flow field. The convection-diffusion equation is 



where p is fluid density, Pi is pressure, and ti is dynamic viscosity. 

For computations on single-crypt/villus units, as shown in Fig. 2 A, 
constant inhibitor concentration was imposed at the outer boundaiy 
and an insulation boundary condition (no flux through membrane) 
was imposed at the crypt/ villus surface. Volume flux at the crypt sur- 
face in the steady state, Jv(z), depends on volume flux in the absence 
of inhibitor, Jv°(z), and local inhibitor concentration, Ci(ro,z), 



JvW= rW/[l + Ci(ro,z)/X,], 



(4) 



where Kd is the equilibrium dissociation constant for inhibitor-target 
(CFTR) binding. 

To model a membrane-permeable inhibitor in which inhibitor is 
uansported from the circulation into the crypt/viUus lumen, the 
transepithelial flux of inhibitor, Jb(z), is described by 



Pinh - [Cb-C,(r,z)]= Jb°(z) . [l-Q(ro,z)/Cb], (5) 



where Pi„b describes inhibitor permeability across the intestinal epi- 
thelium, and Cb is inhibitor concentration in blood. A convective 
flux boundary condition was imposed at the outer boundary and at 
the crypt/villus surface. 

hi multi-crypt computations, the same flux boundary condition is 
imposed at the crypt surface. The parabolic velocity profile at the 
lumen inlet was calculated from mean lumen velocity as 



U(x,y,z) = 2U_„- [l-(s/So)2], 



(6) 



where So is the distance from the lumen center to the top of villi, and 
s is the distance from the center of the lumen. Symmetric boundary 
conditions are imposed at the two side walls and a slip boundary 
condition at the top surface. Details on boundary conditions are pro- 
vided in the Appendix (Figs. Al and A2, and Table 2) . 

The computation time to obtain the steady-state solution was 
~10 min for single-crypt computations and 4—24 h for multi-crypt 
computations, as performed on an HP Z600 workstation (Xeon 
E5645 CPU and 32G RAM; Intel) . The time step was automatically 
computed in COMSOL Multiphysics from mesh size and property 
value variation. Computation validation studies are provided in the 
supplemental text. 

Inhibition of net fluid secretion was computed as the ratio 
of total crypt-villus fluid secretion in the presence versus absence 
of inhibitor. 



% inhibition = 100 • (1 "Ziv/ZJ,")- 



0) 



For modeling the kinetics of inhibitor washout, the fraction of 
bovmd inhibitor, fb (0 to I, fb = C-,/ (Ci + IQ) in the steady state), is 
described by the differential equation, 



dfb(z)/dt = ki- [l-fb(z)] ■ C, 



(8) 



at 



convection diffusion 



(1) 



which contains convection and diffusion terms. The velocity field, Vi, 
in the convection term was computed from the Navier-Stokes equa- 
tion for an incompressible fluid and the continuity equation, 



pf^ + p(v.-v) 



-VPi +riV^Vi 



(2) 



where k_i is inhibitor dissociation rate constant, and ki is inhibi- 
tor bimolecular association rate constant, svibject to the condi- 
tion, = k_i/ki. 

The percent inhibition of fluid secretion along the length of 
the intestine was deduced from the computational results done 
for small segments, in which inhibitor concentration in each 
segment is diluted progressively becavise of fluid secretion. Con- 
servation of inhibitor molecules (Eq. 9) and of fluid volume 
(Eq. 10) requires 

(9) 



Jin et al. 263 



U„,.-AL=Ui„-A,+2:j,- i^d^/^), (10) 

where Cin and Com are inhibitor concentration at the inlet and 
the outlet in each segment, Uin and Uom are mean lumen veloc- 
ity, Al is luminal cross-sectional area, dL is lumen diameter, w is 
segment length, and EJ, (equal to EJv° [1 — 0.01 ■ % inhibition]) 
is total single-segment secretion rate. Percent inhibition at each 
segment was specified by an empirical fit of computed results at 
the segment. 

Total percent inhibition (over the length of the intestine) is 
the ratio of integrated secreted fluid without versus with inhibi- 
tor. Total secreted fluid in the absence of inhibitor was computed 
by summation of flux from each segment, 

(11) 

Secreted fluid (no inhibitor) = X! J ° ' ("-"^l / w) ■ (Lint / '^^)> 

where IJ, ° is initial single-segment secretion rate, dL is lumen 
diameter, w is segment length, Lin, intestinal length, and Ax is 
segment length. Total secreted fluid in the presence of inhibi- 
tor was computed from the difference in luminal flow in the 
final segment without versus with inhibitor, so that percent in- 
hibition becomes 



total % inhibition = ^ 

100 • [l-(Uon,'"'-Uin')AL/l2:j;-('tdL/w)-(Ln„/Ax)|], 



where Uout'* is mean flow velocity in the final segment, and Uin' is 
mean inlet flow velocity at the first segment, where Uqu,^ is calcu- 
lated iteratively applying conservation conditions (Eqs. 9 and 10). 

Online supplemental material 

The online supplemental material includes additional informa- 
tion on the mathematical modeling, including the validations of 
computation time and the number of mesh elements, and details 
of the secretion ratio (],°) calculation. It is available at http:// 
www.jgp.org/cgi/content/full/jgp.20 1 2 10885/DCl . 



RESULTS 

Single crypt-villus computations for a membrane- 
impermeant, surface-targeted inhibitor 
We first analyzed convection-diffusion for a single crypt- 
villus unit to explore the main features of the model. 
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Figure 2. Convective inhibitor washout in 
single crypts. Inhibitor convection-diffu- 
sion was solved for a cylindrical ci-ypt-villus 
unit exposed at its outer boundary (intesti- 
nal lumen) to a constant concentration of 
membrane-impermeant inhibitor, Co/J^. 

(A) Cross section of a cylindrical crypt- 
villus unit showing relative J,° at indicated 
regions (left) . Two-dimensional compu- 
tation volume showing mesh elements 
(middle) . Equivalent three-dimensional 
crypt-villus unit showing downward in- 
hibitor diffusion and upward convection. 

(B) Steady-state profiles of drug concen- 
tration (left), fluid flow velocity (middle), 
and pressure (right) for single-crypt diffu- 
sion-convection for indicated Co/ Ad, with 

= 4 X 10-2 pL/cmVs. (C) Steady-state 
profiles for indicated Jv°, with Co/Ki = L 
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Fig. 2 A (left) shows a cross section of the crypt-vilhis 
structure in mid-jejunum (where the majority of fluid 
secretion occurs in cholera) , which consists of a 350-]am- 
long villus shown aligned with a 150-pm-long crypt. 
Based on reported data on CFTR distribution (Jakab 
et al., 2011), relative XJi ^roiZ) in the absence of inhibitor 
was taken as 1.0 (lower 50% of crypt), 0.6 (upper 50% 
of crypt), 0.28 (lower 50% of villus), and 0.15 (upper 
50% of villus) . Fig. 2 A (middle) shows mesh element in 
the two-dimensional computation region, which, by sym- 
metry, solves the equivalent three-dimensional crypt- 
villus problem (Fig. 2 A, right) . 

For the two-dimensional computations, the intervil- 
lous region was modeled as a simple cylinder approxi- 
mating the actual three-dimensional geometry, which is 
modeled below. We note that the further simplification 
of our model to one dimension is not possible because 
of the need to impose unrealistic boundary conditions, 
which necessitated the inclusion of a "buffer zone" at 
the crypt-villus outlet in the single-crypt computations. 
As fluid leaves the crypt-villus unit it takes with it solute 
and so produces a complex two-dimensional profile ex- 
tending from the oudet that can be quite large for high 
flow velocities as found in diarrhea. 

Fig. 2 B shows steady-state profiles of drug concentra- 
tion (left) and fluid flow velocity (middle) for different 
(dimensionless) Co/K^, in which Co/ IQ is fixed at the 
outer computational region. The increase in pressure is 
negligible, as shown in the right panel of Fig. 2 B for 
Co/K^ = 100. Convective effects are observed in the pro- 
files, with marked axial concentration and velocity gra- 
dients. At fixed fluid secretion rate (in the absence of 
inhibitor) of Jv° = 4 x 10"^ pL/cm'/s, increasing Co/ 
to 100 allowed inhibitor access deep in crypts (Fig. 2 B, 
left). Increased inhibitor concentration reduced fluid 
velocity. Fig. 2 C shows the steady-state drug concentra- 
tion and fluid flow velocity profiles at fixed Co/ -Kd of 1 
for different Jv°, and the pressure profile for Jv° = 10"^ 
pL/ cm^/ s. Increasing Jv° reduced inhibitor concentra- 
tion in crypt-villus units because of convection, with a 
very small pressure increase of <2 Pa (compared with 
atmospheric pressure of 10" Pa) . 

Results from computations as in Fig. 2 are summa- 
rized in Fig. 3. Fig. 3 A shows the percent inhibition of 
secreted fluid for the crypt-villus geometry of the mid- 
jejunum, as deduced from volume flux exiting the crypt 
without inhibitor versus with inhibitor. The percent in- 
hibition of fluid secretion increased from 0 to 100% 
with increasing Co/ K^, as expected, although the appar- 
ent IC50 for inhibition was substantially increased with 
increasing Jv°. For low, near static flow (]° = 10"^ pL/ 
cm^/s), Co/-Kd of 1 produced ^^50% inhibition of fluid 
secretion, as expected. For Jv" of 7 x 10 ' pL/cm^/s, 
which is typical of cholera, inhibition of net fluid secre- 
tion by 50% required Co/ K4 of ^200, and clinically rel- 
evant inhibition by 90% required Co/ -Kd of ^750. Fig. 3 B 



shows similar computations for the colon, which con- 
sists of long, narrow crypts (length 430 pm, diameter 
16 pm) without villi. A similar dependence of percent 
inhibition of fluid secretion on Co/ IQ and Jv" was seen, 
in which 50% inhibition of fluid secretion required 
Co/Ka of --40 for J,° = 7 X 10"^ pL/cmVs, typical of 
cholera in colon. 

Single crypt/villus computations for a membrane- 
permeable, surface-targeted inhibitor 
Computations were also done for a membrane-permeable, 
surface-targeted CFTR inhibitor in which the inhibitor en- 
ters the crypt-vUlus lumen from the blood compartment 
rather than the intestinal lumen. We assume constant in- 
hibitor concentration in blood, with entry into the crypt- 
villus lumen determined by the (blood-to-lumen) inhibitor 
concentration gradient and enterocyte permeability. 

A mid-jejunum 
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Figure 3. Convective inhibitor washout requires a very high con- 
centration of a membrane-impermeant inhibitor in the intesti- 
nal kimen for antisecretory efficacy. (A) Computations done for 
human mid-jejunal anatomy as in Fig. 2. Percent inhibition of net 
secreted fluid as a function of Co/IQ for indicated J,". Jv° of ^7 x 
10"^ pL/cmVs is typical in cholera. (B) Computations done for 
human colonic anatomy with long, narrow crypts without villi. 
Percent inhibition of net secreted fluid as a function of Co/ IQ 
for indicated J,°. J,° of ^7 x 10"^ jiL/cm^/s in colon is typical 
in cholera. 
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The hypothesis is that high blood inhibitor concentra- 
tion and enterocyte permeability are required for anti- 
secretory efficacy of a membrane-permeant, externally 
targeted inhibitor. Our reasoning, as diagrammed in 
Fig. 4 A, is that inhibitor molecules crossing the entero- 
cyte barrier are convected away from the luminal sur- 
face, impairing their ability to inhibit CFTR. 

Fig. 4 B shows drug concentration profiles for differ- 
ent drug permeability coefficients for transepithelial 
transport from blood into the crypt-villus lumen, Fi^h 
(left) , and for different fluid secretion rates, Jv°. Con- 
vective washout of luminal inhibitor is seen as reduced 
C;/ Cb with increasing J,,". Fig. 4 C summarizes the per- 
cent inhibition of net secreted fluid as a function of Pinh 
for different Jv° at fixed Cb = 10 pM and Kd=lO ]iM. At 
a lowJv° of 10^^ pL/cm^/s, half-maximal inhibition re- 
quired Pinh of ^10"^ cm/s, whereas at high Jv° of 7 x 
10^^ pL/cmVs, typical of cholera, 50% of maximal in- 
hibition required an ^100-fold increased Pjnb of ^10^^ 
cm/s. Maximal inhibition was 50% for Cb = K^, as ex- 
pected. Similar computations are shown in Fig. 4 D for 
Cb = 1 pM and -Kd = 0.1 pM, where maximal inhibition is 
^100%. At lowJv° of 10"^ pL/cm^/s, half-maximal inhi- 
bition required Pinh of '^lO^® cm/s, whereas at higher 
Jv° of 7 X 10"^ pL/cm"/s, half-maximal inhibition re- 
quired Pinh of '^5 X 10"^ cm/s. These computations 
show that convective washout effects can occur for a 
membrane-permeable, surface-targeted inhibitor. 



Multi-crypt computations with three-dimensional 
intestinal geometry 

Computations were also done for the three-dimensional 
intestine as diagrammed in Fig. 1 B. Our approach was 
to first analyze percent inhibition versus Co/ -Kd (as done 
in Fig. 3 A for single crypt-villus units) for a short seg- 
ment of intestine (mid-jejunum), and then extend the 
computation to model many meters of intestine. Fig. 5 
A shows an en face view of the crypt-villus geometry in 
this model, as deduced from electron microscopy data 
(Loehry and Creamer, 1969; Marsh and Swift, 1969). Fin- 
ger-shaped villi of 100-250-pm diameter extend out of the 
plane, with inter-crypt distance of ^100 pm and a crypt to 
villus ratio of ^3. We modeled crypts and villi as cylinders 
with semispherical outer surfaces. Villus length weis taken 
as 350 pm, with an outer diameter of 120 pm and inter- 
villtis spacing of 1 50 pm. Crypt length weis taken eis 1 50 pm, 
with an inner diameter of 20 pm and inter-crypt distance 
of 100 pm. The villi to crypt ratio was 3, giving a crypt 
density of 13,333 cm ~, similar to the measured value of 
12,500 cm ' (Trbojevic-Stankovic et al., 2010). Fig. 6 A 
(middle) shows a small segment of 1-mm length and 
150-pm width used as the computational volume (chosen 
to allow a crypt to villus ratio of 3) . The computational 
volume is composed of two layers, mucosa and lumen, of 
different viscosities. Fig. 5 A (right) shows 80,000-150,000 
mesh elements in the computation volume. 

Fig. 5 B shows the steady-state profiles of drug con- 
centration. Similar to the profiles in Fig. 2 B, increasing 
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Figure 4. Model of antidiarrheal 
efficacy of a membrane-permeant 
inhibitor. Computations done for 
human mid-jejunal crypt-vilhis 
anatomy. (A) Schematic show- 
ing inhibitor permeation across 
enterocytes from blood into the 
crypt-villus lumen, which de- 
pends on inhibitor transepithe- 
lial permeability coefficient, Pi„h, 
and blood concentration Cb- 
(B) Steady-state inhibitor concen- 
tration profiles for indicated Pinh 
forJ,° = 2 X 10"^ pL/cmVs (left), 
and for indicated J,° for Pi^i, = 
10"^ cm/s (right), both for Cb = 
10 pM and = 10 pM. (C) Per- 
cent inhibition of net secreted 
fluid as a function of Pinh for in- 
dicated J," at Cb = 10 pM and = 
10 pM. (D) Same as in C, with Cb = 
1 pMand/5i = 0.1 pM. 
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Jv° produced greater convection, reducing inhibitor 
access deep in crypts. Fig. 5 C shows percent net fluid 
inhibition versus Co/-Kd for different ]°, anafogous to 
the single crypt-villus analysis in Fig. 3 C. Co/ was '^1 
for low Jv°, as expected, increasing substantially for 
higher Jv°. These computations for three-dimensional 
crypt-villus geometry support the robustness of the con- 
clusion that greatly increased C0/K4 is needed for effec- 
tive inhibition at high fluid secretion rates. 

The short segment of intestinal surface with anatomi- 
cally accurate three-dimensional geometry was extended 
to model the human mid-jejunum with 5-mm length and 
5-cm diameter (Fig. 5 D, top) . Assuming a well-mixed in- 
testinal lumen beyond the level of villus tips, the steady- 
state profile of fluid inhibition along the intestine was 
computed as described under Model computations in 
Materials and methods. Fig. 5 D (bottom) shows percent 
inhibition along the intestine for different initial lumen 
axial flow velocities, Umean, representing different axial 
starting points along the length of the intestine. The per- 
cent inhibition decreaises along the length of intestine as 
inhibitor is diluted by secreted fluid from crypts and villi. 



Fig. 5 E shows deduced percent inhibition versus Co/ 
plots for different Umean for Jv° of 7 x 10"^ \iL/ cvc?/ s, typi- 
cal of cholera. Because of inhibitor dilution effects, the 
percent inhibition versus Co/-Kd profiles is right-shifted, 
indicating the requirement of high Co/-Kd (>1,000) for 
efficient reduction in intestinal fluid secretion. Even 
greater Co/Ki would be predicted without the assump- 
tion of a well-mixed intestinal lumen. 

Inhibitor washout kinetics limits antidiarrheal efficacy 

The computations described above for an orally admin- 
istered (membrane-impermeable) inhibitor were done 
for the steady state, with constant inhibitor concentra- 
tion entering the intestinal lumen (at the boundary of 
the computation region) . However, because of intermit- 
tent oral administration, inhibitor concentration varies 
over time, with litde inhibitor present in the intestinal 
lumen much of the time. To model inhibitor washout, 
inhibitor concentration in the lumen was reduced at 
time 0 from Co/ -Kd to 0 (step-function) . Fig. 6 A shows 
the kinetics of inhibitor washout for a single crypt-villus 
unit for dissociation rate constant, k_i = 10^^ min^^. 
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Figure 5. Inhibitor washovit in a 
tliree-dimensional multi-crypt model 
of the intestine. (A) Crypt-villus 
geometry in the three-dimensional 
model showing a matrix of spatially 
distinct narrow crypts and relatively 
wide villi. (B) Drug concentration 
profiles for steady-state solution 
of multi-crypt diffusion-convection 
for indicated J," for CJfQi = 10. 

(C) Percent inhibition of net se- 
creted fluid as a function of Co/ IQ 
for indicated J,°. J," of --7 x IQ-^ 
pL/ cm^/ s is typical in cholera. 

(D) Percent inhibition along the 
intestine for indicated lumen axial 
mean velocity, Umean- Segmental per- 
cent inhibition in C was described by 
an empirical regression to tire equa- 
tion: % inhibition = a ■ exp[p ■ logio 
(Ci„/Ad)], with oi = 3.523 and 3 = 
1.112 for J,° = 7 X 10-2 jiL/cmVs, 
typical of cholera. (E) Percent in- 
hibition for 5-m-long intestinal seg- 
ment as a function of Co/ Ki for 
indicated Umean- Results for a short 
(1-mm) intestinal segment are shown 
for comparison. 
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Washout occurs over time as CFTR inhibitor dissocia- 
tion releases inhibitor molecules that are convected out 
of the crypt. In addition to crypt-villus geometry, the 
determinants of washout kinetics include k_i, fluid se- 
cretion rate, and Co/Ki- Fig. 6 B shows the kinetics of 
loss of inhibition of fluid secretion for mid-jejunum, 
with parameters chosen to give '^90% inhibition of 
fluid secretion before washout. For a k_i of 10 ' min^' 
(well below that expected for small-molecule inhibitors 
with micromolar K^) and a high Co/ of 10, washout 
occurs in tens of minutes. Very low k_i is required to 
gready slow inhibitor washout. For the three-dimen- 
sional crypt-villus geometry, predicted washout is even 
more rapid, occurring over a few minutes for a k_i of 
lO^^min^i (Fig. 6C). 



Qualitatively, the effects of convection and diffusion 
on washout can be derived from a simple dimensional 
analysis of the convection and diffusion terms in Eq. 1 , 
giving a single dimensionless parameter, the Peclet 
number (Pe), 



washout effect(Pe) = 



convection term V ■ VC( 



L . 



diffusion term DV 



ViL , L^J ° 
» 4 — — 
D dD 



where Vi is flow velocity, L is length of crypt or villus, D 
is diffusion coefficient, Jv" is initial volume flux, and d 
is crypt/ villus diameter. The washout effect is pro- 
portional to crypt or villus length and volume flux, and 



drug concentration 



DISCUSSION 

CFTR offers a unique target for antisecretory therapy of 
enterotoxin-mediated diarrheas because it is expressed 
on the luminal membrane of enterocytes and is hence 
exposed to the aqueous-phase contents of the intestinal 
lumen. The luminal exposure of CFTR allows the pos- 
sibility of developing inhibitors that are targeted to its 
external surface that could act after oral administration 
and without systemic absorption. High throughput 
screening of chemically diverse small-molecule libraries 
identified glycine hydrazides and related analogues that 
blocked the external CFTR pore direcdy (Muanprasat 
et al., 2004). The glycine hydrazides were engineered 
into nonabsorbable macromolecular conjugates with 
high affinity binding to CFTR and slow dissociation 
(Sonawane et al., 2007, 2008). However, because of the 
high volume of fluid that can be secreted in diarrhea, 
inhibitor washout from the surface of the intestine is a 
concern. Convection-diffusion, as fluid is secreted out 
of the long, narrow intestinal crypts, is predicted to 
reduce inhibitor concentration at the luminal sur- 
face of the crypt epithelial cells. Washout effects can 
also be important for systemically absorbed inhibitors 
that target an extracellular site on crypt epithelial 
cells, as inhibitor dissociation from the surface and 
subsequent convection out of crypts can reduce surface 
inhibitor concentration. Although orally administered 
nonabsorbable drugs are in use (Charmot, 2012), such 
as cholesterol-binding resins and antimicrobials, drug 
targeting to the external site of a luminal receptor on 
the intestinal epithelium is a new concept in drug deliv- 
ery. Because there was no precedent for existing drugs 
using this targeting strategy, or data to evaluate the 
magnitude or kinetics of washout effects, we modeled 
the CFTR-inhibitor system in the intestine using data 
applicable to secretory diarrhea, although our model 
would apply in general to drugs that target luminal re- 
ceptors on enterocytes. 
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Figure 6. Rapid washout of a luminally delivered membrane- 
impermeant inhibitor. After steady-state inhibition, inhibitor 
concentration in the intestinal lumen was reduced from Co/ to 0. 
(A) Drug concentrations profiles for kinetic (presteady-state) 
solution of single-crypt diffusion-convection for Co/IQi = 10, Jv° = 
10"^ pL/cmVs, and k_i = 10"' min"'. (B) Kinetics of inhibition 
of fluid secretion after inhibitor washout. Percent inhibition of 
net fluid secretion as a function of time after washout for Co/ IQl = 
10, J,° = 10"' pL/cmVs, and indicated k_i. (C) Washout kinetics 
for the three-dimensional multi-crypt model for same parameters 
as in A. 
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TABLE 1 

Material properties and baseline parameters 



Material property Baseline parameter 

p Fluid density = 10^ (kg/m') 

Dl Diffusion coefficient in the intestinal lumen = 2 x lO"'* (m^/s) 

Dm Diffusion coefficient in the mucosa (crypt/villus region) = 2 x 10"^"^ (m^/s) 

T7l Fluid dynamic viscosity in the lumen = 10"^ (Pa • s) 

T^M Fluid dynamic viscosity in the mucosa = 10"^ (Pa • s) 

Ly Villus length = 350 (yim) 

Lc Crypt length - 150 (j-mi) 

dy Villus diameter - 120 ((im) 

dc Crypt diameter = 20 (|.im) 

dL Intestinal lumen diameter = 5 (cm) 



inversely proportional to crypt or villus radius and 
diffusion coefficient. Substantial washout effect occurs 
for small crypt/villi radius and long crypt/villi length. 
The geometric term (L^/d) for villi (L = 350 pm, d = 
120 pm from effective diameter) is similar to that for 
crypts (L = 150 pm, d = 20 pm) , which gives a Pe of ^14 
and 16, respectively, for the severe diarrhea case {\" = 7 
X 10 " pL/cm^/s). Therefore, both crypts and villi are 
important determinants of washout, although the wash- 
out effect in crypts is greater because of higher fluid 
secretion rate. In any case, quantitative analysis required 
full computational analysis because of nonlinear effects 
and the complex geometry of crypt-villus units. 

The main conclusion of the modeling is that the anti- 
secretory efficacy of a surface-targeted channel inhib- 
itor requires both (a) very high luminal inhibitor 
concentration compared with its K4, and (b) sustained 
very high luminal inhibitor concentration or slow in- 
hibitor dissociation compared with the interval between 
oral dosing. For a crypt secretion rate of 7 x 10^^ pL/ 
cm^/s, typical of cholera, the luminal concentration of 
a membrane-impermeant inhibitor should be >200 x 
to reduce net fluid secretion by 50% and 800 x to 
reduce fluid secretion by 90% (Figs. 3 A and 5 C) . Even 
higher inhibitor concentration is needed if the inhibi- 
tor diffusion coefficient is lower than that used in the 
model here, which may be the case because of mucus 
and other viscous substances in the aqueous-phase 
fluid. Further, even if adequate inhibitor concentration 
in the intestinal lumen is achieved for antisecretory ef- 
ficacy, a high concentration must be sustained, as deliv- 
ery of inhibitor-free fluid at the luminal surface results 
in rapid loss of inhibitor efficacy unless inhibitor disso- 
ciation from its luminal target is very slow. For a typical 
small molecule with low micromolar K4, dissociation oc- 
curs over a few seconds or less, producing significant 
loss of antisecretory efficacy over minutes (Fig. 6) . For 
oral administration of an inhibitor every few hours, the 
lumen is exposed to cyclical changes in inhibitor con- 
centration, with little inhibitor present most of the 
time. Sustained release or other formulations, although 



potentially helpful, are unlikely to produce sustained 
therapeutic inhibitor concentrations when the fluid se- 
cretion rate is high as in cholera and traveler's diarrhea. 

Of the available externally targeted CFTR inhibitors, 
GlyH-101 (Sonawane et al., 2006) and small-molecule 
analogues such as MalH (Sonawane et al., 2007) and 
iOWH032 (de Hostos et al., 2011) have IC50 in the 
range of 2-8 pM and dissociate from CFTR within sec- 
onds. We conclude from the modeling here that such 
compounds would not fulfill criteria (a) or (b) . MalH- 
lectin conjugates have higher affinity (-Kj of ^50 nM) 
and slow Wcishout over several hours or longer (Sonawane 
et al., 2007), as the lectin moiety is bound tighdy and 
trapped in the dense glycocalyx lining the intestinal sur- 
face. However, although MalH-lectin conjugates can thus 
fulfill criteria (a) and (b) , lectin conjugates may not be 
practical for antisecretory therapy in developing coun- 
tries because the lectin moiety significantly increases 
the drug cost and reduces its stability during storage. 

The antidiarrheal efficacy of a surface-targeted, mem- 
brane-permeable inhibitor that is systemically absorbed 
and enters the intestinal lumen from the blood side was 
also considered. Fluid convection can cause washout of 
inhibitor molecules that enter the crypt lumen from the 
blood. Typical intestinal transepithelial permeability 
coefficients for absorbable drugs, based on Caco-2 data, 
are in the range of 10^® to 10 cm/s, with values up to 
7 X 10 cm/s for the most highly permeable drugs 
(Usansky and Sinko, 2005). The computations in Fig. 4 
indicate that antidiarrheal efficacy is limited by inhibi- 
tor convective washout at high drug Cb for typical 
drug permeability (10^^ ^ 10 ' cm/s). A more potent 
drug (IQ of ^0.1 Cb) predicted good antidiarrheal effi- 
cacy. Therefore, provided that therapeutic drug con- 
centration in the blood can be maintained, and that 
drug Xd is sufficiendy low and drug permeability suffi- 
ciendy high, a surface- targe ted, membrane-permeable 
inhibitor can be efficacious in secretory diarrhea. 

Although our model was formulated to closely reca- 
pitulate the in vivo three-dimensional geometry of a 
high density of radially oriented crypt-villus units lining 
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the intestinal lumen, several features of the intestine 
were not included in the model. The model does not 
consider differences in the detailed geometry and trans- 
port properties of the different regions of the small and 
large intestine, nor does it include fluid absorption or 
mechanical peristalsis. Although these additional com- 
plexities were not modeled explicitly, the major con- 
clusions of our model are robust and should not be 
influenced by these fine details. The majority of se- 
creted fluid comes from the jejunum, as was modeled 
here, and fluid secretion in secretory diarrheas such as 
cholera gready exceeds absorption. Peristaltic effects 
would be predicted to have litde time-averaged effect 
on inhibitor diffusion in the unstirred crypt luminal 
compartment. 

Notwithstanding these potential limitations of our 
model, we conclude that washout in secretory diarrhea 
is a major factor for nonabsorbable drugs that target 
surface receptors or transporters on enterocytes in in- 
testinal crypts. Therapeutic efficacy requires sustained 
high drug concentration in the intestinal lumen. For 
the inhibitor-CFTR system, surface-targeted small 
molecules with micromolar are unlikely to be effica- 
cious for therapy of secretory diarrhea, whereas mac- 
romolecular conjugates with nanomolar IQ and very 
slow dissociation might be effective. We conclude, there- 
fore, that although surface-targeted glycine hydrazide 
and MalH originally appeared to be attractive candi- 
dates for antisecretory therapy, absorbable CFTR in- 
hibitors, such as thiazolidinones (Thiagarajah et al., 
2004) and PPQ/BPO compounds (Tradtrantip et al., 
2009; Snyder etal., 2011), are better development can- 
didates for CFTR inhibition therapy of enterotoxin- 
mediated secretory diarrheas. 

APPENDIX 

Boundary conditions 

As described in Materials and methods, the three- 
dimensional multi-crypt computation involves solution of 
the diffusion-convection equation for inhibitor con- 
centration and the (incompressible fluid) Navier-Stokes 



equation for the fluid field. Fig. Al and Table 2 provide 
details on boundary conditions. 

Lumen inlet boundary condition 

The circular lumen inlet has parabolic velocity pro- 
file with mean velocity, Umean- Because the intestinal 
lumen (^5-cm diameter) is much bigger than crypt 
dimensions (20-]am diameter and 150-|im length), 
the computation was simplified to greatly reduce 
computation time. A small segment near the mucosa 
(crypt and villus), with computation volume of 1,250 
pm (crypt length of 150 pm, villus length of 350 pm) , 
was used (Fig. 1 C) . The velocity profile at the lumen 
inlet was taken as part of the parabolic profile as 
shown in Fig. A2. The lumen inlet boundary condi- 
tion (the magnitude of mean velocity) had negligible 
effect on percent inhibition (not depicted) because 
the viscosity in the mucosa is > 10-fold greater than that 
in the intestinal lumen, so that the flow field produced 
by secretion is dominant in the crypt/villus mucosa. 
A constant inhibitor concentration boundary condi- 
tion was used for the inhibitor concentration at the 
lumen inlet. 

Left and right surface boundary condition 
As shown in Fig. 5 A, the small segment containing 
crypts and villi unit is repeated symmetrically around 
the lumen, producing laminar flow in the lumen (vis- 
cous forces dominant over inertia) , R^ = Umem'i-L/ (t) < 
0.1, where dL is lumen diameter and v is kinematic vis- 
cosity. Thus, symmetric boundary conditions were im- 
posed at the left and the right surface boundaries. 

Top surface boundary condition 

A slip boundary condition was used for the Navier-Stokes 
equation, and an insulation (no flux) boundary condi- 
tion Wcis used for the diffusion-convection equation. 

Crypt and villus surface boundary condition 

As described in Materials and methods, the rate of 
fluid secretion at the crypt and villus surface was deter- 
mined by initial secretion rate and inhibitor concentra- 
tion. For a membrane-impermeant inhibitor insulation 



TABLE 2 

Boundary conditions for Navier-Stokes and diffusion— convection equations 



Navier-Stokes equation Diffusion-convection equation 



(1) Inlet BC 
(lumen flow velocity) 


U = 2U„ea„ ■ [1-(S/S„)^] 

V,= -U-n 


(1) Inlet BC 
(constant C^) 


C 


= c„ 






(2) Symmetric BC 


n-Vi = 0, t • [-pi + p(VVi + (Wi)^] = 0 


(2) Symmetric BC 


n 


■ n = 0, N 


= -D 


■ VC + CVi 


(3) Slip BC 


n-Vi = 0, t • [-pi + p(VVi + (Wi)''] = 0 


(3) Insulation BC 


n ' 


■ ;z = 0, N 


= -D 


■ VC + CVi 


(4) Flux BC 


V, = -Jvn,J,,=J.°(l -C|/(Q + ;^,)) 


(4) Insulation BC 


n ' 


■ n = 0, N 


= -D 


■ VC + cv, 


(5) Outlet BC 
(no viscotis stress) 


[Vi(VVi+ (VV,)^]n = 0,p = p„ 


(5) Outlet BC 
(convective flux) 


n- 


(-DVC) 


= 0 





BC, boundary condition; n, surface normal vector; I, unit vector; V, gradient operator. Other variables are defined in the main text. 
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A Navier-stokes equation 



B diffuion-convection equation 




Figure A1 . Schematic of boundary conditions for the three-dimensional multi-crypt computation. (A) Boundary conditions for the 
Navier-Stokes equation. (B) Boundary conditions for the diffusion-convection equation. 



(no flux) , boundary conditions were imposed at the crypt 
and villus surfaces. 
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Lumen outlet boundary condition 

Because the lumen outiet is open, a no-viscous stress and 
pressure boundary condition Wcis used for the Navier- 
Stokes equation, and a convective flux boundary condi- 
tion Wcis used for the diffusion-convection equation. 

This work was supported by National Institutes of Health grants 
DK72517, EB00415, HL73856, DK35124, and EY13574. 

Edward N. Pughjr. served as editor. 



lumen 



computation volume 




U(x,y,z,t) = 2U„ea„ |l-(s/s„)^l 



Figure A2. Schematic of inlet boundary condition, where s is dis- 
tance from the center of the lumen, So is distance from the center 
of the lumen to the top of the villus, and Umean is the mean inhibi- 
tor velocity in the lumen. 
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